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ABSTRACT 

One of the most dramatic signatures of the reionization era may be the enormous ionized bubbles around lu- 
minous quasars (with radii reaching ^ 40 comoving Mpc), which may survive as "fossil" ionized regions long 
after their source shuts off. Here we study how the inhomogeneous intergalactic medium (IGM) evolves inside 
such fossils. The average recombination rate declines rapidly with time: the densest pockets recombine rapidly, 
leaving low-density regions that recombine much more slowly. Furthermore, the brief quasar episode signif- 
icantly increases the mean free path inside the fossil bubbles. As a result, even a weak ionizing background 
generated by galaxies inside the fossil can maintain it in a relatively highly and uniformly ionized state. For 
example, galaxies that would ionize 20-30% of hydrogen in a random patch of the IGM can maintain 80-90% 
ionization inside the fossil, for a duration much longer than the average recombination time in the IGM. Quasar 
fossils at z < 10 can thus retain their identity for nearly a Hubble time, and will appear "gray," distinct from 
both the average IGM (which has a "swiss-cheese" ionization topology and a lower mean ionized fraction), and 
from bubbles around active quasars (which are fully ionized). More distant fossils, at z > 10 have a weaker 
galaxy-generated ionizing background and a higher gas density. They can break up and attain a swiss-cheese 
topology similar to the rest of the IGM, but still with a smaller contrast between the ionized bubbles and the 
partially neutral regions separating them. Analogous Helll-fossils should exist around the epoch of Hell/Helll 
reionization at z 3. Rapid recombinations inside the Helll-fossils will be more common, because many of 
them will have no Hell-ionizing background; but the time lag before another quasar appears is typically small, 
so < 50% of the gas is able to recombine, even for fossils that form well before reionization is complete. 
Our model of inhomogeneous recombination also applies to "double reionization" models and shows that a 
non-monotonic reionization history is even more unlikely than previously thought. 
Subject headings: cosmology: theory - intergalactic medium 
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1. INTRODUCTION 

The reionization of hydrogen and helium throughout the in- 
tergalactic medium (IGM) are landmark events in the early 
history of structure formation. As such, they (and particu- 
larly hydrogen reionization) have received a great deal of at- 
tention - both observationally and theoretically - in the past 
several years. At present, t he observational evidence leaves 
significant ambiguities (see iFan et al.1l2006al for a recent re- 
view). The electron scattering optical depth of r ^ 0.09 
measured through cosmic microwave background polariza- 
tion anisotropics implies that hydrogen reio nization began at 
z > 10, albe it with a large uncertainty (Page et al.l 120071: 
iDunklev et al.l l2008l) . On the other hand, Lya absorption 
spectra of quasars at z ^ 6 show some evidence for a rapid 
transitio n in the globally-averaged neutral fraction, xm (e.g., 
iFan et al. 20Q6b ). Howev e r the L ya absorption is so saturated 
in the ,Gunn & PetersonI (Il965h trough (with optical depth 
tqp > lO^xm) that constraints derived dir ectly from the lack 
of flux in that spectral region are weak (" White et all 120031: 
iFan et al. 2006b). Studying the properties of the ionized zones 
and the flux transmitted at the corresponding wavelengths 
near the Lya line has yielded much tighter lower limits on 
the neutral fraction along a few sight lines, using Ly a and 
Ly/3 transmission statistics (Mesinger & Haiman 2004, '2007^ 
and the inferred ionized zone sizes (Wyithe & Loeb 2004a; 
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iMesinger & Haim ar? '2004), although the interpretations of 
these resu lts is stil l subject to poor statistics and possible bi- 
ases (e.g, iLidz et all 2006: Becker et all 120071: iMaselh et all 
l2007t iBolton & Haehnelt 2007b.a) . Helium reionization is 
equally controversial: althoug h Hell Lya forest spectra sug- 
gest r apid evolution at z '--^ 2.9 dHeap et al.l2000l:ISmette et all 
l2002h . secondary signatures (such as the IGM thermal his- 
tory and radiation background) suggest change at slightly 
higher reds hifts - or in some cases no e volution at all (e.g., 
Schave et al . 2000; Theuns et al. 2002; F aucher-Gig uere et al] 
2007^, 'Songaila 1998; Kim et al. 2002t lAguirre et al.. i2004l: 
Spngaila 2005; Bernardi et al. 2003). 

One of the most interesting aspect of the reionization pro- 
cess is its inhomogeneity: different regions of the Universe 
can be reionized at very different times, depending on their 
relative proximity to the ionizing sources. During reioniza- 
tion, the resulting patterns of ionized and neutral gas can then 
be a rich source of information about t he IGM, the ionizing 
sources, and the ways they interact (Furlanetto et al.l l2004t 
iFurlanetto &'Ohll2007bl) . In essence, large ionized zones pin- 
point either large clusters of sources or exceptionally bright 
ones. 

The most dramatic examples are luminous quasars, which 
can ionize enormous regions, spanning several tens of (co- 
moving) Mpc, in either hydrogen or helium. These sources 
are (probably) rather short-liv ed, with expected lifetimes ^ 
10''-10^ yr (see, e.g.. iMartinil 12004 for a review). When the 
quasar shuts down, the ionizing background within the region 
drops precipitously: the gala xies that surr ound the quasar host 
may be relatively numerous dAlvarez & A bel 2007; iLidz et all 
l2007t iGeil & Wyiflieir2007h but generically cannot compete 
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with such a strong ionizing source. As such, the IGM will 
begin recombining as soon as it shuts off. 

However, the recombination time of IGM gas (at the mean 
density) is comparable to the Hubble time at z ~ 8 and 
z ^ 3 for hydrogen and helium, respectively. Thus dur- 
ing either recombination era one expects the ionized re- 
gion to survive for a relatively long time after its source 
vanishes, becoming a "fossil" bubble. As some of the 
largest coherent features in the Universe at these high red- 
shifts, these bubbles (both during their active and fossil 
phases) are obvious objects of interest. In particular, dur- 
ing hydrogen reionization they are favored targets of first- 
generation 21 cm surveys, which can only make images on 
the largest scale s - and so are limited to search for structure s 
of these sizes ( Wvi the & Loe b '2004b'; IWvithe et alj l2005h . 
The characteristics of the se regions can be used to study the 
prope r ties of the quasar ([Wvithe et al.l 120051: IZaro ubi & SiU? 
2005"; 'Thomas & Zaroubil l2008t iKramer & Haiman. ,2007i: 
Sethi & Haiman 2008) as well as its re l ation to the sur- 
rounding s ources ( Alvarez & Abell l2007t iLidz et all 120071: 
iGeil & Wvi the 2007). Although active bright quasars are 
quite rare at high redshifts, the long lifetime of fossils will 
make large features much easier to find - increasing their 
abundance by approximately the ratio of the recombination 
time to the source lifetime. During helium reionization, the 
distribution of these fossils will affect the thermal histories 
of IGM gas elements and hence their Lya forest properties 
dGleser et al.ll2005l:lFurlanetto & Ohll2007al) . 

Despite their promise, there has been relatively little atten- 
tion paid to the physics of the fossils, and in particular to how 
they actually recombine. The most naive picture, in which the 
IGM is assumed to be uniform at the mean density, is obvi- 
ously too simple because, even at these high redshifts, the fil- 
amentary structure of the cosmic web is growing - and it has 
long been recognized that the accompanying clumping is cru- 
cial for understanding the reionization process. In that con- 
text, the most common approximation for this inhomogeneity 
is the clumping factor C = (ri?^ / {rif' , which is the ratio 
of the recombination rate (volume-averaged over all ionized 
regions) to its naive expectation in a uniform medium."* 

This simple approach can be adequate so long as the spa- 
tial averaging to derive the clumping factor excludes neutral 
regions, which are of course not recombining, while simulta- 
neously accounting for the large-scale structure in the IGM 
dMiralda-Escude et alJlIOOOllFurlanetto & Ohll2005h . In that 
case, so long as the division of the IGM into neutral and 
ionized regions is not changing rapidly, a constant clumping 
factor approach is useful because the densest ionized regions 
(which account for most of the recombinations) are relatively 
stable. 

However, with fossil bubbles this assumption is manifestly 
false: after such a precipitous drop in the ionizing back- 
ground, the entire region falls out of ionization equilibrium. 
The dense gas will initially be highly ionized by the quasar, so 
the clumping factor will be large. But these regions recombine 
quickly and become neutral, making the total recombination 
rate (and effective clumping factor) fall, eventually flattening 
out at a level appropriate to the smaller ionizing background 
from the remaining galaxies. 

In practice, it is often approximated by averaging over t he entire volume , 
excluding only regions bound to dark matter halos (e.g., Illiev et al] 120061 : 
iMellema et al. 2006), because weigh ting by the ionized fraction intr oduces 
extra complexity (but see iGnedin & Ostrikeill997l : IKohler et aU2007r) . 



Similar physics should describe any global "recombination 
era", during which the global average ionized fraction is de- 
creasing with cosmic time. This may occur in "double reion- 
ization" scenarios, in which there is a sharp transition in the 
mode of star formation so that the IGM is reionized twice, 
with substantial recom binations in between (.Cen 2003b.al : 
IWvithe & Loebll2b03allbh . although in pra ctice such histories 
are difficul t to arran ge self-consis tentlv jHaiman & Holdeii 
2003; Furla netto & L oeb 2005; Iliev et alJl2007h . In this case, 
the entire IGM would be "fossilized", with the second genera- 
tion of stars unable to maintain the same level of ionization as 
the previous one. Such models have only been studied using 
a simple clumping prescription, and they rely on a relatively 
large value for C to speed up recorabina tions during the inter- 
mediate phase (Furlanetto & Loeb 2005). But the same argu- 
ments apply as for fossils: a large clumping factor will only 
be appropriate until the dense gas recombines. Afterward, the 
effective clumping factor will decrease, and the global recom- 
bination rate will slow down. 

The goal of this paper is to describe the physics of this re- 
combination process in some detail and thus to evaluate bet- 
ter the structure of fossil ionized regions. We will aim to 
answer several basic questions. Does the dumpiness of the 
IGM allow most of the gas to recombine, or only the densest 
regions? Can the residual ionizing background from galax- 
ies maintain high ionization, even if they could not provide 
the initial ionization of the gas on their own? How does the 
ionization topology evolve inside the fossils? Does the in- 
homogeneous IGM allow regions far from galaxies to shield 
themselves and recombine faster than those near the galaxies 
- so that the fossil resembles the swiss-cheese topology of the 
rest of the IGM? Or, does the ionization remain more uniform 
inside fossils than outside? And, finally, can quasar fossils be 
distinguished, either individually or statistically, from large, 
galaxy-generated ionized bubbles, and from active quasars, 
through the topology of their ionized gas and their abundance? 

This paper is organized as follows. In §|2]and[3] we describe 
our models for computing the recombination histories of gas 
parcels and the density distribution of the IGM. Then, in fQ- 
|6l we apply this model to fossils produced during hydrogen 
reionization, studying, respectively, a fossil with no ionizing 
sources, a fossil with a uniform ionizing background flux, and 
an alternative, "global" fossil model, incorporating the inho- 
mogeneous density field of the IGM, but explicitly conserving 
photon number In § |2l we apply these models to fossils dur- 
ing helium reionization. Finally, in § [8] we discuss some of 
the observational prospects to detect these fossils, and we of- 
fer our conclusions in § |9l 

In our numerical calculations, we assume a cosmology 
with = 0.26, f^A = 0.74, f]f, = 0.044, Hq = 
100ft. km s-i Mpc-i (with h = 0.74), n = 0.95, and 
( 78 = 0-8, consistent with the most rec ent measurements 
jDunklev et alJl2008l:lKomatsu et al.ll2008l) . Unless oflierwise 
specified, we use comoving units for all distances. 

2. THE RECOMBINATION HISTORY 

Our first task is to compute the recombination history of 
an individual gas par cel. We will use a siniplified v ersion of 
iHui &Gnedi"iy(II997h : see lFurlanetto & OhI (l2007al) for more 
details. The full solution requires us to follow the abundances 
of three independent species (including hydrogen, helium, 
and electrons), the gas temperature (which affects the recom- 
bination and collisional ionization rates), and the density of 
the parcel. 
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Consider a gas element of fractional overdensity S and tem- 
perature T illuminated by an ionizing background of the form 



- 



1 I^HI < < t'HcII 
/ i^HoII < 1^, 



(1) 



where Jhi.-21 is the angle-averaged specific intensity of the 
background, in units of 10^^^ cm^^ s^^ Hz~^ sr~^, z/hi and 
i-'Hcii are the frequencies corresponding to the HI and Hell 
ionization edges, and / is a constant. We will assume for 
simplicity that / = before helium is reionized and f — I 
afterward. In our fiducial models, we adopt a = 1.5 (typi- 
cal of quasar spectra). After the quasar that initially ionized 
a region turns off, the ionizing background will be provided 
by stellar sources, and the spectrum will soften considerably. 
However, the shape only changes the photoheating rate, which 
does not substantially affect our calculations. (It does not af- 
fect the Helll fraction because we set / = before helium 
reionization anyway.) 

We define the number density of species i to be ri; = 
(1 + 5)XiPb/rap, where pf, is the mean proper mass density 
of baryons. Note that the Xi are the species number fractions, 
not the neutral fractions. 

The thermal evolution of the element is determined by 



dt 3{i + d)dt 



dQ 



dt 



SkBTih dt ' 
(2) 

where d/dt is the Lagrangian derivative. The first term on the 
right hand side describes the Hubble expansion, the second 
describes adiabatic cooling or heating from structure forma- 
tion, the third accounts for the change of internal energy per 
particle from changing the total particle density, and in the last 
term dQ / dt is the net heat gain or loss per unit volume from 
radiation processes (see below). 

The second term requires an expression for the growth of 
the nonlinear density field. iHui & GnedinI (11997 *) used the 
Zel'dovich approximation, along with an analytic estimate for 
the distribution of the strain tensor, which has been shown to 
provide an accurate analytic approximation to the density evo- 
lution in the quasilinear regime. However, to maintain com- 
putational simplicity, we will assume that the gas elements 
evolve following the spherical collapse (or expansion) model. 
We map the linear densiti es to nonlinear over densities via the 
following fitting formula (iMo & WhitellI996h . 



Sl = D{z)Sl = S,- 



1.35 



1.12431 



0.78785 



(l + <5)2/3 ' (1 + (5)0-58661 ' 

(3) 

where D{z) is the linear growth function and 6c ~ 1.69 is 
the threshold for virialization in the spherical collapse model. 
Under this assumption, we can write 



dS_ _ . dP dS 

dt " " 'dfds2' 



(4) 



In the present context, this approximation affects the temper- 
ature distribution (and hence recombination coefficient) sub- 
stantially only at densities far from the mean, so it is not a 
significant concern for us. 

The fourth term includes a number of radiative heating and 
cooling processes. The most important heating mechanism is 
photoionization itself; each species i contributes a term 



dQ 
dt 



(5) 



where ai is the photoionization cross section for species i, Vi 
is its ionization threshold, and h denotes Planck's constant 
(to di fferentiate it from the Hubble constant). We use the 
fits of IVerner et al.l (|l996) for the photoionization cross sec- 
tions. The other relevant mechanisms cool the gas, and in- 
clude Compton cooling off the CMB (which dominates dur- 
ing the era of hydrogen reionization), recombinations (radia- 
tive and, for He II, dielectronic), collisional ionization, col- 
lisiona l line excitation, and free-free emission. We use the 
fits of iHui & GnedinI (Il997h for all of these processes ex- 
cept Compton coo ling (for which we use the exact form in 
ISeaeer et al.llI999|) and free -free emission (for which we use 
the fit in lTheuns et alJl998l) . We have verified t hat our results 
are un changed if we use the fits pr esented in iTheuns et al.l 
(!1998) (an updated form of those in ICen|[T992l) for all these 
mechanisms. 

Finally, we have an ionization balance equation for each 
species; e.g., for hydrogen. 



dX, 



HI 



dX 



HII 



Tm+af{T)XcXmi[nb{l+5)], 



(6) 



dt dt 
where as is the case-B recombination coefficient and 

r,: = / dv " o-,: (7) 



is the ionization rate of species i. We will use T12 — 
r/(10~^^ s~^) for con venience. Note that, unlike in 
iFurlanetto &~Ohl (l2007ah . we do not assume photoionization 
equilibrium for these calculations. 

3. THE DENSITY DISTRIBUTION 

Equations (|2]), (|4|i, and (|6|l allow us to follow the evolution 
of a single gas parcel. We wish, however, to consider the ag- 
gregate evolution of large regions of the IGM; we therefore 
need the distribution of gas density inside the IGM. Unfortu- 
nately, there is no well-tested model at high redshifts. In the 
following we will use the density distribution Pv(A) (where 
the distribution is over vo l ume a nd A = 1 + 5) recommended 
bv lMiralda-Escude etalj (l200d henceforth MHROO), which 
fits cosmological simulations at z = 2-4 quite well: 



Py(A) dA = AoA^'^exp 



(A-2/3 - C^f 



2(2(5o/3)2 



dA. (8) 



Intuitively, the underlying Gaussian density fluctuations are 
modified through nonlinear void growth and a power law tail 
at large A. MHROO argued that the form could be extrapo- 
lated to higher redshifts in the following way. First, (5o essen- 
tially represents the variance of density fluctuations smoothed 
on the Jeans scale for an ionized medium (appropriate for our 
purposes, where we will only consider gas that has been pre- 
ionized); thus 5^ cx (1 + z)^^. The power-law exponent (3 
determines the behavior at large densities; we set (3 — 2.5 
for z > 6. The remaining constants (^0 and Co) can be 
set by demanding proper mass and volume normalization. 
Note, however, that this distribution comes from a simula- 
tion with somewhat different cosmological parameters than 
the currently preferred values and so is probably not accurate 
enough for detailed comparisons. 

MHROO also offer a prescription for determining A;, the 
mean free path of ionizing photons, in a highly-ionized uni- 
verse. They assume that the gas density field has two phases: 
low-density gas is highly ionized and transparent to ionizing 



4 



photons, while high density regions (with A > A^) are neu- 
tral because of self-shielding. We will examine the accuracy 
of this assumption la ter, but it will be quite useful for most 
of our work (see also iFurlanetto & Ohl l2005). In this picture, 
the mean free path equals the mean distance between the self- 
shielded clumps along a random line of sight, which is ap- 
proximately 

A, = Ao[l~i^y(AO]-'/^ (9) 
Here Fy ( A j ) is the fraction of volume with A < A^ and Aq is 
a (redshift-dependent) normalization factor. Formally, this ex- 
pression is valid only if the number density and shape (though 
not total cross section) of absorbers is independent of A^ . This 
is obviously not true in detail for the cosmic web. However, 
MHROO found that equation (|9]l provides a good fit to numer- 
ical simulations at z = 2-4 with XqH{z) — 60 km s^^ (in 
physical units). We will extrapolate the same prescription to 
higher redshifts. Equation (|9]l was derived in the highly ion- 
ized limit, ignoring photoelectric absorption by low-column 
density systems, and therefore generically overestimates Ai 
at the ionization edge. The correction ca n be up to a fac- 
tor of ~ 2 even in highly ionized cases ( IFurlanetto & OhI 
l2005h . Some of this decrease can be compensated by the 
longer mean free paths of higher-energy photons, but to be 
conservative, we set \qH{z) — 30kms^^ whenever we re- 
fer to the MHROO mean free path in our calculations. We will 
discuss the accuracy of this mean free path estimate in more 
detail below. 

The MHROO prescription was developed for the post- 
overlap universe, when dense regions are self-shielded and the 
mean-free path is controlled by Lyman-limit systems. In the 
pre-overlap universe, assuming that all regions with A < A,; 
are ionized implicitly assumes an "outside-in" reionization 
topology, in apparent contradiction with expectation s from 
semi-analytic models and s imulations (e.g., Furlane tto et al.l 
l2004t iMcOuinn et al.l |2007|) that reionization is an "inside- 
out" process due to source bias. The MHROO prescription is 
likely applicable during reionization once the mean free path 
is sufficiently large that a typical photon samples a fair (and 
relatively un biased) fraction o f the density distribution dur- 
ing its flight (IZhang et al.ll2007h .^ The MHROO prescription is 
particularly well suited to quasar fossil bubbles. For one, the 
initial mean free path within the quasar bubble (tens of Mpc) 
is comparable to that of the post-overlap universe, when Fy is 
close to unity; one is considering a sufficiently lar ge volume 
that it is a fair sample of the density distribution. iLidz et al] 
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( l2007h find that the bias is modest when averaged over the 
quasar bubble volume. Moreover, dense regions recombine 
first, leaving underdense regions ionized — ^precisely the as- 
sumptions made in the MHROO model, provided the ionizing 
background is relatively uniform (see f|6T| for more discus- 
sion of the latter point). 

Figure [T] shows A; as a function of mass-averaged ion- 
ized fraction Xi = F^i^i) for several different redshifts 
(z = 2, 3, 4, 6, 8, 10, and 15, from top to bottom). Note that 
Xi decreases as redshift increases because the IGM gets less 
and less clumpy, so the spacing between gas above a given 
density threshold decreases. At higher redshifts, Xi w 0.5 
corresponds to A w 1, so these regions are typically sepa- 
rated by a fraction of a Mpc. 

It is important to note that the MHROO mean free path is 
much larger than one would expect in a uniformly ionized 

' Note that the details of the small-scale density distribution within and 
around halos are generally absorbed into the escape fraction. 




Fig. 1 . — Comoving mean free path of photons at the hydrogen ionization 
edge as a function of the ionized fraction Xi in the MHROO model. The curves 
assume z = 2, 3, 4, 6, 8, 10, and 15, from top to bottom. 



medium. The optical depth in a uniform medium, measured 
across one fiducial MHROO mean free path, can be expressed 
as 



T\-^xmnHKo-m 



(10) 



' 260 a;Hi 



3 Mpc 



1 + z 
9 



where xhi is the mean neutral fraction of the gas outside of 
self-shielded regions, hn is the mean density of hydrogen nu- 
clei, and (Thi is the photoionization cross section. For the 
latter, we used (Thi = (t = 2 x lO^^'^ cm^, the valiie averaged 
over a typical stellar spectrum (Miralda-Esc udell2003i) . Fig- 
ure [T] shows that under the MHROO geometry, the mean free 
path reaches ^^Mpc scales when xm 0.1 — 0.5. In com- 
parison, in a uniform IGM, reaching this large mean free path 
would require .xhi ^ 10^"^. The reason is that the MHROO 
procedure calculates the mean free path in an entirely dif- 
ferent way: it assumes that the bulk of the IGM is highly- 
ionized (with negligible absorption) except for regions above 
the self-shielding threshold (which are entirely neutral). In 
other words, in the MHROO geometry, xi (which is a mass av- 
erage) refers to the fraction of gas outside these self-shielded 
regions, not to the volume-averaged ionized fraction. We will 
check for the consistency of this picture later and see some 
examples in which it becomes untenable. 

4. FOSSIL BUBBLES WITH ZERO IONIZING BACKGROUND 

We begin by considering the simple case of an ionized bub- 
ble whose generating source turns off instantaneously, with 
no sources remaining or turning on. Thus the gas simply re- 
combines, without any further influence from a radiation field. 
In the context of hydrogen reionization, pre-existing galaxies 
in the fossil bubble will always maintain an ionizing back- 
ground. Considering the zero-background case is neverthe- 
less a useful pedagogical exercise; furthermore, in the con- 
text of helium reionization, where the death of the quasar may 
leave no ionizing sources, this case is also physically relevant. 

To produce our ionization histories, we initialize our cal- 
culation at an assumed redshift Zi where the quasar has fully 
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Fig. 2. — Ionization liistories for fossil HII regions that are produced when 
the ionizing source turns off at some initial redshift Zi . The fossils are as- 
sumed to see no ionizing flux, so the gas is allowed to recombine. The three 
sets of curves assume Xi = 1 at Zi = 8, 10, and 15, from left to right. 
Within each set, the solid and dashed curves show the mass- and volume- 
averaged ionized fractions, Xi^^n and „, respectively, using the MHROO fit 
for the density distribution. The dotted curves show the histories for elements 
at the mean density (with A = 1). 



ionized a region. We ignore the residual neutral gas at this 
stage, including any self-shielded, high-density clumps that 
may remain neutral even in the presence of a quasar. We then 
compute Xi{A,z) at a final redshift z. The mass-averaged 
ionized fraction within the bubble, Xi,m, is 

^vm = /iGM/ dAx,(A,z)AFy(A). (11) 

We set Amax = 50 so as to exclude virialized gas from the 
calculation (since its ionization properties are irrelevant to 
the IGM properties). The prefactor /jo'^ is the fraction of 
mass in gas with A < A^ax and normalizes the result to 
give the ionized fraction of just the IGM gas. The volume- 
averaged ionized fraction, Xi^i, is computed similarly, except 
with APy(A) — ^ P\/(A) inside the integrand and /igm the 
volume fraction of gas with A < A,„ax- 

Figure |2] shows the evolution of the ionized fraction in a 
fossil bubble with = 0. We consider three initialization 
redshifts, Zi ~ 8, 10, and 15. For each case, we show Xi^m, 
Xi^y, and Xi{A — 1) by the solid, dashed, and dotted curves 
(respectively). As one would expect from the short recom- 
bination time at these high redshifts, all of these cases are 
able to recombine significantly. Also as one would expect, 
Xi^m < Xi^y, because underdense regions (which recombine 
slowly) of course fill more than their share of the volume. 
However, the two are not so far apart, and in particular the 
mass-averaged value is always quite close to the expectation 
for gas elements at the mean density. 

Many models (including both semi-analytic and numeri- 
cal varieties) compute reionization histories using a so-called 
clumping factor, C = (n\^ / {rie)'^ to estimate the total re- 
combination rate. The dotted curve in Figure |3] shows this 
factor in the MHROO model as a function of redshift, assum- 
ing that all gas with A < 50 is ionized. The effective value is 
always larger than unity, but not by an extremely large value; 
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Fig. 3. — The effective clumping factors for the recombining fossils in 
Fig.|2] For reference, the dotted curve shows C for the MHROO IGM model, 
assuming that all gas at A < 50 is fully ionized (with gas above this threshold 
residing inside collapsed halos and excised from the IGM). 



it is significantly smaller than often assumed because we have 
excluded virialized systems and also assumed that the gas is 
Jeans-smoothed on the scale appropriate for ionized gas. It 
increases toward lower redshifts as structure formation con- 
tinues and amplifies density contrasts. 

Unfortunately, blindly applying this clumping simplifica- 
tion only works when all the gas is in ionization equilib- 
rium. In reality, this is never a particularly good assump- 
tion. For example, during reionization, the dense interiors of 
self-shielded clumps remain neutral, hardly ever see ionizing 
photons, and do not contribute to the ionized volume. Ion- 
izing photons need not battle recombinations in this gas, but 
only the slower recombinations in the more rarefied outskirts 
of such regions. The dotted curve essentially deals with this 
problem by excluding dense gas from the clumpi ng calcula- 
tion (as is typical in nume rical simulations; e.g., Illiev et al.l 
120061; iMcOuinn et al.ll2007i) . In reality, one should allow the 
maximum density threshold A^ to evolve to reflect the actual 
ionizing background (and neutral fraction) within the bubble. 

The problem is somewhat more insidious for fossil bubbles. 
The solid, long-dashed, and short-dashed curves in Figure [3] 
show the effective clumping factors for fossil bubbles, again 
with initial ionization at zi = 8, 10, and 15. Here we define 
the clumping factor Ccfr by comparing the overall recombina- 
tion rate to that for a gas element at the mean density. After 
an initial spike caused by rapid cooling of the dense gas (re- 
member a (X T^^ ''),^ Ccff quickly decreases to a value close 
to unity. 

This behavior occurs because there is (by definition) zero 
ionizing background inside fossils, so no gas elements can be 
in ionization equilibrium. Although dense elements recom- 
bine quickly, they cannot help their less-dense neighbors to 
do so; once a region is mostly neutral, it breaks off from the 

* This spike is probably not physical, because the gas is actually heated 
over the lifetime of the quasar and expands and cools continuously, rather 
than being set to an arbitrary initial temperature; see Shapiro et al. 1 2004); 
IMesinger et alj <2006D for discussions of the thermal behavior of photoheated 
gas in simulations. 
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Fig. 4. — Neutral fraction as a function of density inside fossil HII regions. 
(a): 2i = 15. The curves show z = 14-6, from bottom to top, with Az = 2; 
these have ini.m = 0.32, 0.67, 0.82, 0.89, and 0.93. (b): Zi = 10. 
The curves show z = 9-6, from bottom to top, with Az = 1; these have 
iCHl,m = 0.68, 0.5, 0.62, and 0.70. The filled hexagons mark the nominal 
Ai values, which, in the MHROO model, produce the same average neutral 
fraction. 



recombining gas. Thus the density of an "average" recom- 
bining parcel decreases with cosmic time, and with it the ef- 
fective clumping factor. It eventually levels off at Ccff < 1, 
because the recombination time for gas at the mean density is 
comparable to the Hubble time. This behavior is vital for an- 
alytic models of fossil bubbles - and also for any model with 
a "recombination era", including so-called "double reioniza- 
tion" - and a proper treatment requires a full accounting of the 
IGM density distribution. This also implies that cosmological 
simulations cannot use simple subgrid clumping in any cells 
where the local ionizing background is zero; instead the in- 
stantaneous recombination rate depends on the integrated re- 
combination history of the cell. 

Another way to see this is through the neutral fraction as 
a function of density. Figure |4] shows sequences of a;Hi(A) 
at several different times, with Zi = 15 and Zi = 10 in 
the top and bottom panels, respectively. For the most part, 
low-density gas remains highly-ionized, while dense regions 
quickly recombine. However, the transition between the two 
is rather extended (because the recombination time is propor- 
tional to A^^). Indeed, even deep voids can have a substantial 
neutral fraction long after the quasar turns off. 

4.1. The Mean Free Path 

Next we will consider how the MHROO model for the ion- 
ized gas distribution and mean free path fares in this scenario. 
Qualitatively, one might expect it to perform adequately: af- 
ter all, the distribution goes from relatively highly-ionized gas 
when A ^ 1 to nearly neutral gas at high densities. In 
fact, blindly using the ionized fraction to fix A; (employ- 
ing the physical scenario of MHROO, and assuming that gas 
with A < Ai is actually highly ionized) even roughly repro- 
duces the regime for which Xi ~ 0.5 (see Fig.|4|. However, 
Ai clearly does not mark a sharp separation between highly- 
ionized and nearly neutral gas; rather, the low-density gas is 
also marginally neutral. As a result, the mean free path for 



these fossils will be far smaller than shown in Figure [T] Sim- 
ply assuming a uniform medium (eq. [Tol l yields t ^ 1 across 
a region of size ^ 7 kpc at .thi ^ 0.5 - no larger than a 
typical halo! 

We can do slightly better than this uniform IGM approx- 
imation by dividing the gas into discrete sys tems in a simi - 
lar way to the Lya forest at lower redshifts ( ISchaveir2001h . 
The typical length scale of such a system is the Jeans length, 
Lj ^ G/y/CsP, where Cg is the local sound speed. Then the 
column density of neutral hydrogen (which determines the 
optical depth of each system to ionizing photons) is A'hi ~ 
uhLj cx n]/^T^/^. At the Lyman edge. 
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Comparison to Figure |4] shows that the neutral fractions re- 
quired to produce t > 1 are typically exceeded even in low- 
density gas. Thus, effectively, every discrete absorber - even 
a deep void - is a Lyman-limit system, and the mean free path 
of an ionizing photon is extremely small. The MHROO model 
is obviously not applicable in this scenario. 

5. A NONZERO IONIZING BACKGROUND 

We now include a small, but nonzero, ionizing background 
in our calculations. This background could represent pre- 
existing galaxies, as well as any new galaxies that continue 
to form within the fossil bubble, and/or any low-level resid- 
ual emission from the quasar after it has turned off (see ^16. 31 
below). To set the amplitude during hydrogen reionization, 
we begin by assuming that the emissivity is simply propor- 
tional to the rate at which gas collapses into star-forming ha- 
los, d/coii/di- Here /coii is the collapsed fraction in halos 
above the minimum mass m,ni„ able to form stars. We will 
parameterize this threshold mass in terms of the mass m4 
of halos with virial tempe ratures larger than 10^ K, so that 
atomic cooling is efficient jBarkana & Loebl[2001h 7 We as- 
sume an ionizing efficiency (, such that the globally-averaged 
ionized fraction at the initial redshift Zi is^ 

= Cfcolliz^)■ (B) 

5.1. Gas at the Mean Density 

We can now take two possible approaches. In this sec- 
tion, we will discuss an approach similar to § |4l i.e. track- 
ing the evolution of each gas element independently, simply 
adding a redshift-dependent ionizing background to the evolu- 
tion equations. This approach is imperfect, as we will discuss 
below, but it is useful to provide a feel for the effectiveness of 
an ionizing background in suppressing recombinations. In the 
next section, we will follow an alternative (more global) ap- 
proach, based directly on the emissivities (rather than fluxes). 

We begin by examining the evolution of gas with A = 1. 
The ionization rate can be written as 



- ni,aX 



d(C/o 



dt 



(14) 



(15) 



Note that it is possible for the effective mmin to differ between the 
universe at large and the highly-ionized fossil bub ble, becau se (for exam- 
ple) phot oheating may increase the Jeans mas s Eees 1986; Ef stathioiJll99a 
IThoul & Weinberg., 1996i ; jDiikstra et aU2004[) . 

^ Note that this expression ignores recombinations, but in our calculations 
they can be simply swept into the free parameters described below, so long 
as they are roughly uniform. 
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Fig. 5. — Ionization liistories for gas with A = 1 under a uniform ion- 
izing background. Tlie two sets of curves assume = 1 at = 8 and 
10. Tlie solid curves assume F = as in Fig. |2] The long-dashed, short- 
dashed, short dot-dashed, dotted, and long dash-dotted curves within each set 
(i.e., from top to bottom) take (zj)Ai = 0.6, 6, 60, 600 and 6000 kpc. 
respectively. Note that this calculation breaks down when self-shielding be- 
comes important; see text. 



where ei, is the emissivity per unit frequency and Xi, is the 
(proper) mean free path of a photon with frequency v. Substi- 
tuting values appropriate for our cosmology, we find 



Ti2 = 0.018 s 
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(16) 

where is the comoving mean free path (averaged over fre- 
quency) and we have used equation ( fTST l to eliminate the effi- 
ciency C, in terms of the more physically transparent globally- 
averaged ionized fraction when the bubble is first ionized.^ 

Figure|5]shows ionization histories for gas at the mean den- 
sity, initially ionized at z,; = 8 and 10. The solid curves 
assume F = 0. The long-dashed, short-dashed, short dot- 
dashed, dotted, and long dash-dotted curves within each set 
(from top to bottom) take xl{zi)\i ~ 0.6, 6, 60, 600 and 
6000 kpc, respectively. This combination of parameters sets 
the overall amplitude of the ionizing background in equa- 
tion (O; the background can be large either because one be- 
gins late in reionization (if is large) or because the mean free 
path in the fossil bubble is large. For comparison to results 
below, it is useful to think of these as having {zi) = 0.2, so 
that the bubble is created early in reionization (but not too 
early) and A^ = 0.003-30 Mpc. Comparison to Figure [T] 
shows that the MHROO model could apply to the IGM in 
the latter two cases (3-30 Mpc); the other cases imply very 
small mean free paths and correspond to a regime where the 
MHROO assumption breaks down. 

5.2. The Mean Free Path 

We will next explicitly consider how the MHROO model for 
the mean free path fares in this picture. We must immediately 

' In fact, bright quasars form in the most massive halos and hence in biased 
regions of the IGM, inside of which /coii should be larger than its global 
mean. [Lidz et al. (2007) have shown, however, that although the bias is large 
within a few Mpc of the quasar, it is relatively modest when averaged over 
the entire volume of a large fossil bubble, so we ignore it for simplicity. 



note that Figure|5]shows the evolution only for gas at the mean 
density. However, in the interesting limit in which this gas re- 
mains highly ionized, the fossils approach an asymptotic ion- 
ized fraction relatively quickly: this is the value appropriate 
for photoionization equilibrium, which states 

-3/2 



mi 



(1 - xm) 



1.5 X lO^^A 

fcon{zi) 



Mpc 



l + z 



d/cou/dz 



(17) 



Here we have assumed T = 10** K; in fact, the temperature 
evolves as the universe expands, so this expression is only 
approximately correct. Except for the brief initial adjust- 
ment phase, the equilibrium value works extremely well for 
xl{zi)\i > 60 kpc (so long as the appropriate temperature 
is chosen). Thus we expect xm et A in other gas parcels, 
so long as the ionizing background illuminates every region 
equally. 

Now let us suppose that a quasar appears at z^, ionizes a 
large bubble (with comoving radius > 30 Mpc), and turns 
off. For a concrete example, let us assume that x^ (z^) — 0.2. 
The fossil will immediately begin recombining. At this ini- 
tial stage, the mean free path must be at least as large as the 
quasar bubble (otherwise a single quasar could not have ion- 
ized the entire region). Thus the bottom curve in Figure|5]will 
be most appropriate, and the neutral fraction will initially be 
~ 3 X 10^^. Substituting into equation (fTOl i, the optical depth 
across one mean free path is ^ 0.8. Thus the mean free path 
would remain large, and the entire fossil would clearly remain 
highly ionized. 

Such long mean free paths are quite reas onable at high red- 
shifts, even outside of fossils. For example. iLidz et aP (l2007h 
show that extrapolations of the Lya forest to z > 6 yield mean 
free paths at the ionization edge of ^ 8-15 Mpc. Moreover, 
as described above, the MHROO model predicts Xi > 10 Mpc 
for the highly ionized gas immediately after the quasar turns 
off. 

With the help of equation ( fTTI l, we can also see that denser 
gas will also remain quite ionized: in the example above, 
A = 50 gas would initially have xhi ~ 0.015. However, such 
an approach ignores radiative transfer, which is crucial for our 
purposes (it is instead similar to the way many cosmological 
simulations mimic a uniform radiation field by illuminating 
each and every gas particle identically). The basic point is 
that dense gas elements are not isolated: rather they are part 
of discrete regions, surrounded by other dense elements. The 
accumulated opacity of the outskirts of such objects reduces 
the effective ionizing background in the densest elements - or 
in other words they become self-shielded. When this process 
is included, it is possible to approximate the IGM as two- 
phase, with gas at A < A; highly ionized (and eventually 
settling into ionization equilibrium) and gas at A > A^ self - 
shielded and fully neutral (MHROO, iFurlanetto & Ohll2005h . 
It is these self-shielded elements that set the mean free path 
in the MHROO picture, and it is also these regions that the 
simple-minded approach of this section ignores. 

6. FOSSIL BUBBLES AND INHOMOGENEOUS REIONIZATION 

To rectify the shortcomings of the previous section, we note 
that the principal problem is that a uniform ionizing back- 
ground allows too much ionization. Because each gas parcel 
is treated independently, and because the ionizing background 
is not re-adjusted to account for losses due to recombinations. 
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Fig. 6. — Minimum global ionized fraction required for the background 
ionizing sources to counteract recombinations in fossil bubbles. The thick 
curves use the MHROO clumping factor, while the thin curves assume C^ff = 
1. The solid curves assume rrinjin = m,4, while the dashed curves assume 
'"min = 10m4. The figure suggests that fossils can be kept highly ionized, 
at a level well above the ionized traction in the rest of the IGM. 



photons are not conserved in such a scheme. In more concrete 
terms, the scheme above not only puts photons in the outskirts 
of dense regions but also in their centers; self-shielding de- 
mands that the outskirts consume all the photons, so adding 
more to the center must "over-ionize" the universe. 

We will therefore follow a more global approach to the evo- 
lution of the ionized fraction by enforcing photon conserva- 
tion as our starting point. We first ask what emissivity is 
required to overcome the recombinations within an ionized 
bubble. As in equation ( fT3] l. we parameterize the emissivity 
in terms of the equivalent global ionized fraction outside the 
fossil. 

Figure |6] shows aif the minimum xf required to over- 
come recombinations in a fully-ionized bubble (in this con- 
text, "fully-ionized" implies that all gas elements below — 
50 are ionized). Note that if the emissivity was constant at 
all times, then „j„ would simply equal tgai/^rcc, the ratio 
of the time elapsed since the onset of galaxy formation to the 
instantaneous recombination time. In practice, the emissiv- 
ity increases steeply with cosmic time, which decreases the 
value of j^in- Note further that the recombination rate 
is proportional to the ionized fraction, so the minimum re- 
quirement for a particular bubble should actually be multi- 
phed by its mean ionized fraction at the appropriate redshift. 
We show results for two clumping factors: the upper thick 
set of curves take the MHROO value and are appropriate for 
the period shortly after the fossil forms, while the thin lower 
curves assume C = 1. The latter is appropriate for evaluat- 
ing whether relatively old fossil bubbles that have partially re- 
combined will continue to do so once the dense gas is already 
neutral (see Fig. O. The solid curves assume m„iin = TO4, 
while the dashed curves assume m,„in = 10m4. The lat- 

Unlike in the previous section, here we take the global ionized fraction 
at each redshift, not at the moment of the bubble's creation. Also note that our 
conversion from emissivity to ignores recombinations in the background 
medium, so in reality the required value will be even smaller 



ter, higher mass threshold presents an even weaker require- 
ment on xf , because more massive halos evolve more quickly 
over cosmic time, so that the emissivity that can keep the fos- 
sils ionized corresponds to an even lower background ionized 
fraction. 

We see that, even taking a conservative clumping factor, 
fossil bubbles can only recombine when xf < 0.3-0.4. With 
the MHR clumping value, the minimum increases slowly as 
redshift decreases as the IGM becomes clumpier. If, on the 
other hand, we look at the later stages when C ~ 1 (or roughly 
Xi < 0.75), the requirement decreases to 5^ > 0.1-0.2. In the 
latter case, the minimum decreases as redshift decreases be- 
cause the mean IGM density (and hence recombination rate) 
falls as well. This simple calculation suggests that quasar bub- 
bles will only recombine significantly if they are created well 
before reionization is underway, at < 0.1-0.3. 

To give some concrete scenarios, we track the ionized frac- 
tion within the fossil via 



dxj _ ^ d/coii 
dz dz 



dz 



(18) 



where and Up are the the mean electron and proton den- 
sities at redshift z. We assume (unless stated otherwise) that 
'Timin = ™4 to Calculate /cou- For the effective clumping fac- 
tor, we will generally use the MHROO value Cmhr but also 
compare to a case with C = 1 (where the recombination rate 
is slower). This assumes that the same stellar sources that set 
xf [zi) continue to illuminate the fossil. Of course, if the fos- 
sil persists for long enough, and if quasars are relatively com- 
mon, it is possible that another quasar will form in the same 
region; in that case, the fossil will quickly become completely 
ionized again. 

Figure I2] shows the resulting ionization histories of the fos- 
sil bubbles, in comparison to those of the background, x^ (z), 
for a series of assumed emissivities (for the latter, we can also 
use eq.[T8] but with different initial conditions; we use Cmhr 
for the background, but it makes little difference). In each 
panel, the solid and dotted curves are for fossil bubbles with 
C = Cmhr and C — 1, respectively. The dashed curves show 
xf . In the left panel, Zi = 10, with the three curves in each set 
taking (z;) = 0.4, 0.2, and 0.1, from top to bottom. In the 
right panel, Zi = 15, with the three curves in each set taking 
xf{z,) = 0.02, 0.01, and 0.005, from top to bottom. We now 
discuss each of these cases in turn. 

6.1. A Typical Quasar 

As seen in Figure |7] the requirement that reionization com- 
pletes by 2: = 6 forces x^{zi) > 0.15 for Zi < 10. So, 
for quasars at Zi = 10, the contrast between the background 
and the fossils (or, indeed, between the background and active 
quasar bubbles) is never more than a factor of a few. Unless 
supermassive black holes form long before reionization, this 
case is therefore typical of most fossils. Figure [T] shows that, 
as expected, the ionizing background inside the bubble effi- 
ciently suppresses recombinations - in no case does Xi fall 
below ^ 0.8 for these fossils. This is true regardless of the 
clumping prescription - while the MHROO model allows a 
more rapid initial decline, the two prescriptions converge once 
the dense gas recombines and the clumping factor approaches 
unity. 

Our description is still, however, incomplete, because we 
have not answered one important question: as they recom- 
bine, do the fossils re-develop the "swiss cheese" topology 
typical of reionization, with islands of ionized gas separated 
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Fig. 7. — Top: Comparison of the ionization histories for fossil bubbles and the universe as a whole. In each panel, the solid and dotted curves are for fossil 
bubbles with C = Cmhr and C = 1, respectively. The dashed curves show (with C = Cmhr)- I^^ft- = 10. with the thi'ee curves in each set taking 
(zi) = 0.4, 0.2, and 0.1, from top to bottom. Right: Zi = 15, with the three curves in each set taking (zi) = 0.02, 0.01, and 0.005, from top to bottom. 
Bottom: Mean free paths of ionizing photons in the fossil and in the background universe. The solid curves show from the MHROO model; the dashed curves 
show a fiducial clustering length for galaxies according to {z) (see text for discussion). The parameters are the same as in the top panel, except that we do 
not show the cases with C = 1. 



by a sea of neutral (or in this case partially neutral) hydro- 
gen? Up to this point, we have assumed that "shadowing" 
and other radiative transfer effects are unimportant, so that 
photons are delivered to any region that is not self-shielded 
as needed. Of course, in reality the ionizing sources are clus- 
tered, which modulates the local emissivity. During reioniza- 
tion of the background u niverse, this strongly a ffects the dis- 
tribution of bubble sizes (iFurlanetto et al.ll2004l) . Is the same 
true for recombining fossils? 

To answer this, we will need to compare the mean free path 
within the fossil to the typical scale of fluctuations in the ion- 
izing background generated by the clustering of the ionizing 
sources. First we consider the fossils. Our primary tool for 
computing Ai is the MHROO model. Given an ionized frac- 
tion, this sets a threshold density above which gas is as- 
sumed to be self-shielded, and below which it is highly ion- 
ized. In our case, can be determined by demanding that 
the fraction of mass with A < A^ is equal to Xi, with that, Ai 
follows from equation (|9]l. 

The solid lines in the bottom panels of Figure |7] show this 
prescription; they correspond to the same models as in the up- 
per panels, all using Cmhr for the clumping." Obviously, 
the mean free path falls immediately after the bright quasar 
turns off;'^ it reaches a minimum simultaneously with Xi and 
then increases again. In our models, we assume that all gas 
with A > 50 remains neutral, so the mean free paths even- 
tually match onto Ai corresponding to this density threshold 
(the turnover in the solid curves when Xi w 1). 

The mean free path typically falls to A^ ^ 5-10 Mpc be- 
fore increasing again. If we consider the uniform radiation 

" Note that we use Aq = 30kms^^ here, as described in ijl] This is 
a conservative lower bound on the mean free path, because of high-energy 
photons, so if anything we overestimate the probability that regions break off 
from the uniform background. 

In actuality, the initial mean free path is probably no larger than the fos- 
sil, so the decline may not be quite so severe except for the brightest quasars. 



field models shown in Figure |5] this is a rather large value 
and [for the xf{zi) assumed in this panel] the resulting ion- 
izing background should easily maintain a high level of ion- 
ization in the gas. This is important, because the gas between 
the self-shielded regions must remain optically thin for the 
MHROO model to apply. To examine whether this is the case, 
we show in Figure [U? the neutral fraction as a function of 
density at several stages during this recombination process, 
assuming a uniform ionizing background and neglecting self- 
shielding. The solid, dashed, and dotted curves show .thi (A) 
at z = 7, 8, and 9, respectively (note that the entire universe is 
reionized by z = 6). We take x^ {zi) — 0.2 for all the curves, 
so they correspond to the middle case in Figure [7] We have 
also used the mean free paths and emissivity from this model 
as input in the calculation, so it provides the actual a;Hi for all 
the gas that lies below the self-shielding threshold. The filled 
hexagons show this threshold A^ according to the MHROO 
model (and again taking the recombination history in Fig.]?). 

Obviously, even though the mean free path has declined 
significantly, the low-density gas remains quite highly ion- 
ized throughout the history. At z = 9, 8, and 7, the volume- 
averaged neutral fractions are w (1.9, 0.7, 0.14) x 10^"^ (in- 
cluding only gas with A < A^), respectively. Using equa- 
tion dTol l. these correspond to optical depths across one mean 
free path of tx « (0.9, 0.5, 0.4), from high to low redshift. 

Because these optical depths are somewhat below unity, the 
MHROO model appears to be a reasonable description of the 
IGM inside the fossil: the accumulated opacity of gas between 
self-shielded regions is small. However, the absorption is not 
completely negligible, so there may be some pockets where 
the emissivity is smaller than normal (i.e., voids that are far 
from nearby galaxies), particularly in the earliest phases or 
when x^^{zi) < 0.1, where the optical depths are somewhat 
larger To gauge whether this will occur, we must consider the 
scales over which the emissivity itself fluctuates significantly. 

Before proceeding to do so, we note that Figure [8] also 
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Fig. 8. — Neutral fractions as a function of density at a sequence of red- 
sliifts, assuming tliat the ionizing background and mean free patlis are set by 
the middle of the set of three curves in Fig. [7] . Top: (a) The solid, long 
dashed, short dashed, dash-dotted, and dotted curves take 2 = 8,9, 10, 12, 
and 14, respectively, (b) The solid, dashed, and dotted curves take 2 = 7, 8, 
and 9, respectively. Note that all the xhi curves ignore self-shielding. The 
filled hexagons indicate where the MHROO model predicts that self-shielding 
becomes important. 



lets us check another necessary condition for the MHROO 
model: the self-shielded gas must have had sufficient time to 
recombine and become neutral. This is a particular worry for 
these fossils, because the typical densities of the self-shielded 
"Lyman-limit systems" are not too far above the mean (see the 
filled hexagons in Fig. [8]l. Fortunately, the gas is indeed still 
able to recombine quickly: the recombination time of (fully- 
ionized) gas with A = Ai is, at worst, roughly on-tenth of the 
Hubble time. Thus the self-shielded systems should be able 
to reach Xi ~ 0.1 without any trouble. 

Now we return to the question of whether isolated, low den- 
sity regions will be able to break off and recombine more 
quickly. A necessary, but not sufficient, condition to avoid 
returning to a swiss-cheese morphology and to maintain a 

- _ — 1/3 

uniform ionizing background is \i > dgai = n^^j , where 
rigai is the mean number density of galaxies. If we include all 
halos with T^i, > 10* K, Jgai w 0.6, 0.7, and 0.85 Mpc and 
z = 6, 8, and 10, respectively. In this example, Aj remains 
well above this level, and ionizing photons will, on average, 
be able to reach neighboring galaxies. Thus we would not 
expect the most extreme swiss-cheese topology. 

However, inhomogeneity can actually begin before 
reaches the mean galaxy spacing, because galaxies are highly 
clustered at high redshifts. A convenient way to include 
clustering is to compare A^ to the characteristic size Rc 
of ionized bubbles around clumps of galaxies in the back- 
ground universe. We compute Rc as a function of ionized 
fraction using the analytic model of iFurlanetto et al.) (l2004l) . 
which is accurate to within a factor of ~ 2 d uring most 
of re i onization, and better in the later stages dZahn et all 
12001: iMesinger & Furlanettd l2007h . We show R^ by the 
dashed curves in the bottom panels of Figure |7] (see also 
IFurlanetto et"ani2006ah .'^ 

We cap Rc at the mean free path when = 50 for consistency with 



We see that Xi > Rc throughout the initial decline and 
plateau in Xi. This means that source clustering will be rela- 
tively unimportant while the bubble is recombining; the typ- 
ical scale of variation is smaller than the residual mean free 
path in the bubble, so fluctuations in the ionizing sources will 
be washed out by the long mean free path (as in the post- 
reionization universe). On the other hand, stochastic fluctu- 
ations in the galaxy distribution could still cause some inho- 
mogeneity through a "runaway" effect, in which a region that 
happens to see a lower flux will become more neutral, which 
will help to shield neighboring gas parcels, leading to lower 
fluxes, etc. Such a process would require t\ > 1, which is 
certainly possible in the earliest stages, or if xf{zi) < 0.1. So 
it seems likely that some rare regions far from galaxies will 
begin to recombine. In practice, this runaway effect will prob- 
ably be reduced by the fact that under-illuminated regions are 
already the farthest from the galaxies and so are most likely 
to remain outside the ionized bubbles anyway. 

Once Xi begins to increase, the left panel in Figure|7]shows 
that Xi ^ Rc- (This near equality is coincidental; see the right 
panel for a counter-example.) At this stage, source cluster- 
ing within the bubble will become important. In fact, even 
though the characteristic bubble size in the background uni- 
verse is comparable to or only slightly larger than Xi, the 
"sub-bubbles" within the fossil will actually be somewhat 
larger than in the surrounding universe. 

To understand this, note that Xi reaches a plateau when the 
emissivity is barely enough to overcome the modest recombi- 
nation rate inside the bubbles (we find C w 1 at these times). 
The fossil is then in a "photon-starved" regime where a larger 
emissivity is required in order to ionize additional hydrogen 
atoms. Once Rc ~ Xi, overdense clumps of sources are able 
to do just that over the scale of the mean free path in the back- 
ground universe. But in the fossil, only < 20% of the hydro- 
gen atoms are actually neutral - so a clump of sources in the 
fossil can ionize a volume ~ 5 times larger than an equivalent 
clump in the background universe. The actual sub-clump size 
may even be somewhat larger than this naive prediction, be- 
cause the new, larger sub-bubble can merge with its neighbors 
as well. On the other hand, the quasar is also likely to be in- 
side an overdense region, which would decrease th e expected 
size, a t least in the central, biased parts of the fossil jLidz et al.l 
l2007h . We also note that the actual size of the galaxy sub- 
bubbles can be rather large (several Mpc at z < 8), so that a 
typical fossil may contain only a few dozen such sub-bubbles. 

Meanwhile, if these highly-ionized regions around clumps 
of sources are able to grow because their emissivity is larger 
than average, we also know that underdense voids within 
the fossil will continue recombining even during the phases 
where Xi is either flat or increasing. Thus a contrast be- 
tween the fully ionized regions and their recombining neigh- 
bors will develop, and the fossil will look much like the back- 
ground universe but with larger bubbles (by a factor ^ 1/xi 
in volume) and a significantly smaller contrast between the 
fully ionized regions and the surrounding partially neutral gas. 
Even these shielded regions would not provide much contrast 
because much of the dense gas has already recombined, so 
the effective clumping is modest. We emphasize that, even in 
the absence of any ionizing sources, fossils recombine quite 
slowly in their late stages (see Fig. |2]i, and only a modest 

the fossil treatment. This roughly simulates the appearance of self-shielded 
regions in the background universe; see Furlanetto & Oh (2005J for a more 
complete model. 
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emissivity is required to halt the decrease. (For example, the 
three solid curves span a factor of four in emissivity - much 
larger than we expect from large-scale density fluctuations in 
the IGM - and still the plateau in Xi occurs at a quite large 
value.) 

6.2. An Early Quasar 

Quasars at — 15 will produce fossils with more dra- 
matic contrasts, because reionization at z ^ 6-10 implies 
xf{zi) < 0.02. Thus there is a long time for the bubble to 
remain visible against the mostly neutral background (and, 
eventually, to recombine). However, even in this extreme 
case, the right panels of Figure |7] show that the bubbles never 
become more than about one-half neutral before the back- 
ground catches up to them. One interesting side note is that in 
this high-redshift case, the MHROO clumping factor implies 
more residual neutral gas; as seen in Figure[3] C < 1 except 
in the initial stages, because only underdense gas (where the 
recombination time is long) remains ionized after the initial 
phases. (At z > 14, the MHROO clumping factor does lead to 
slightly faster recombinations, although that is difficult to see 
in the figure.) 

The degree of recombination in these bubbles also suggests 
that the MHROO prescription may break down. Figure[8]con- 
firms this explicitly for the (zi) — 0.01 case (the middle 
curves in the right panels of Figure |7]l. Here we find that 
Ai w 1 for most of the evolution, which leads to mean free 
paths ^ 1 Mpc. (At this extremely high redshift, the IGM 
is still fairly uniform, so there is little variation of the mean 
free path with the ionized fraction; see Fig. [T] Thus falls 
quickly to ~ 1 Mpc and then levels off for a long period, in- 
cluding the plateau and initial increase in Xi.) Until the end 
stages of reionization, the ionized fraction in the A < 
gas (nominally optically thin in the MHROO picture) is a few 
percent at z > 10, so t\ 13, 4.5, 1.5, 0.7, and 0.1 at 
z ~ 14, 12, 10, 9, and 8, respectively. Thus, even assuming 
uniform illumination, the IGM recombines sufficiently to ren- 
der it opaque to ionizing photons on Mpc scales. The MHROO 
picture breaks down, and the mean free path will become a 
fraction of an Mpc. 

Of course, luminous sources will be able to maintain high 
ionization in their own environs, while more distant gas will 
continue to recombine. The mean spacing between galaxies 
at z 15 (10) is 2.1 (0.85) Mpc, respectively. Comparing to 
Figure|2l falls below this level, at least until < 0.5, even 
using the overestimate from the MHROO picture. Thus the 
early stages of the recombination process must be highly in- 
homogeneous: regions surrounding rare galaxies will remain 
highly ionized, but other regions will recombine much like 
the cases shown in Figure |2] As a result, the ionized fraction 
inside of these fossils will be smaller than shown in Figure|7] 
because much of the gas will be recombining without inter- 
ference while the ionizing sources waste photons repeatedly 
ionizing dense gas in their environs. This could lead to much 
higher contrast between the ionized bubbles and partially ion- 
ized surroundings than is ever possible for a fossil surround- 
ing a = 10 quasar. 

Once the clustering of the ionizing sources becomes signif- 
icant, Rc quickly becomes larger than A; (again, even with 
the overestimate from the MHROO picture). This means that 
source clustering will induce stronger variations in the fos- 
sil's ionization structure that begin at a smaller x^ (compared 
to the Zi = 10 fossil). Again, the ionization pattern inside 
the fossil will have a larger characteristic scale than that in the 



background universe, because the surrounding gas has less 
neutral hydrogen, but the contrast between the fully ionized 
bubbles and the rest of the IGM will be somewhat smaller 
than in the background universe. 

Overall, we find that only in extreme cases like this one - 
where the quasar shuts off while the background universe is 
only a few percent ionized - do we expect a well-developed 
swiss-cheese topology inside the fossil. Attaining such a state 
is difficult for two reasons. First, the low-density gas requires 
significant time to recombine; because > 50% of the mass 
has A < 1, the quasar must appear early: at least several 
recombination times before the completion of reionization. 
Furthermore, a swiss-cheese topology also requires that the 
mean free path within the fossil falls below Rc (or, even bet- 
ter, the mean galaxy spacing) so that the clustering pattern 
of the galaxies can induce fluctuations in the emissivity and 
hence ionized bubble pattern. From Fig. |7] we see that A,; 
grows more steeply with Xi than does Rc, only if significant 
recombination takes place does Ai fall below Rc- 

6.3. Residual Emission from the Quasar 

Thus far, we have assumed a "lightbulb" model of quasar 
evolution, in which quasars only turn on for some short pe- 
riod of time, and then shut off, leaving only the remaining 
galaxies to supply an ionizing background to suppress recom- 
binations. However, recent models motivated by hydrody- 
namical simulations of galaxy mergers - in which black holes 
spend the majority of the time well below the peak luminos- 
ity of the associated quasar, but nevertheless radiating - pro- 
vide better fits to multi-wavelength quasar lu minosity func- 
tions and a host of other empirical constraints (iHopkins et al.l 
2005a,b). If the quasar continues shining at some lower lumi- 
nosity, it may suffice to prevent the bubble from recombining, 
even without the aid of neighboring galaxies. If the quasar 
shines at ihi for time thi and then reverts to L\o for an ex- 
tended period, the requirement that the bubble does not re- 
combine during the latter period is (Lio/Lhi) > (ihi/ircc) ^ 

0.03 (thi/10V)C'([l + 2]/10)^ where tree = l/[CaArie] is 
the recombination time with in the bubble. This ac cords well 
with expectations from the (IHopkins et al.ll2005al) model, in 
which quasars shine on cosmological timescales with Edding- 
ton ratios between I ^ 0.01 and Z ~ 1 at peak. A fit to their 
simulations which also agrees with fits to luminosity functions 
is: 

dl^=^Q-P(-^/^Q)' ^19) 

where the fitted parameters Xq = aL-^poak, = 

il^°^(ipcak/[lOi°L0])"^, and aL = 0.20, 4^°^ = 1.37 x 
10^ yr and ut — —0.11. This light curve allows us to cal- 
culate the radiative output as a function of luminosity.'"* For 
instance, although most of the energy is radiated (and hence 
most of the bubble growth occurs) while the quasar is at peak 
luminosity, a Lpoak ~ 10^° quasar spends ^ 1 Gyr with 
L > O.OSipcak, sufficient to prevent the bubble from recom- 
bining until the universe is reionized. In this model, lower 
mass black holes radiate an even large fraction of energy 
at low Edding t on ratio s. Of course, the applicability of the 
IHopkins et al.l ( l2005al) model to high redshifts is uncertain. 

S trictly speaking, in the simulations motivating the IHop kins et alj 
)2005al) fitting formula, a somewhat larger fraction of the low-Eddington ra- 
tio emission is produced before, rather than after the bright, near-Eddington 
phase. 



12 




Fig. 9. — Left panel: Ionization histories for fossil Helll regions with zero ionizing background. The four sets of curves that decrease from right to left assume 
SShcIII = 1 at 2i = 4, 5, 6, and 8. Within each set, the solid and dashed curves show 2'Hcin.m ShcIII,u. respectively. The dotted curves show the histories 
for elements with A = 1. The dot-dashed curve (increasing from the right) shows a sample global reionization history. Right panel: Same as Figure [3] but for 
helium reionization. 



but it is clear that low-luminosity quasar remnants could con- 
tribute very significantly to the ionizing background in 'fossil' 
bubbles. Hence, our arguments considering only the galactic 
contribution provide a strict lower bound on the ionization 
state. 

7. HELIUM REIONIZATION 

7.1. Zero Ionizing Background 

Figure |9] shows corresponding recombination histories for 
fossil Helll bubbles produced during quasar reionization. We 
consider four initialization redshifts here, Zi — 4, 5, 6, and 
8, from left to right, and assume zero ionizing background. 
Qualitatively, the results are similar to those for HI: the mass- 
averaged ionized fraction is typically smaller than the ionized 
fraction for gas at the mean density, but not by a large amount. 
However, once most of the gas has recombined, XHdii.m be- 
comes slightly larger than the ionized fraction of gas with 
A = 1: beyond this point, most of the neutral gas lies in un- 
derdense voids, and the effective clumping factor is substan- 
tially less than unity. The curves are also somewhat steeper 
than for hydrogen because the recombination time for helium 
is several times smaller and because the Universe is clumpier 
at lower redshifts. 

Figure |9] also shows that the clumping factor behaves sim- 
ilarly to fossils during hydrogen reionization (although here 
there is no initial spike, because the temperatures are assumed 
to be hotter at the initial instant of reionization). The asymp- 
totic Ceff are actually smaller here, mostly because the re- 
combination rate is slightly faster relative to the Hubble time 
and so only gas below the mean density remains substantially 
ionized. 

As before, it is useful to compare xhoIII within bubbles to 
the background universe. The dot-dashed curve in the left 
panel o f Figure |9] shows such a history (see Furlanetto & OhI 
2007bl for deta i ls). It uses the quasar luminosity function from 
Hopkins et al.l (|2007), assumes that quasar s have sp ectra typi- 
cal of observed quasars at lower redshifts (IVanden Berk et aTl 
l200ltlTelfer eTal]l2002h and ignores IGM recombinations (so 



is no more than a rough guide). We have reduced the to- 
tal luminosity of each quasar by a factor of ~ 2.5 to force 
reionization to complete at z = 3, consistent with a n um- 
ber of observational probes (see lFurlanetto & Ohl l2007b for a 
summary). Interestingly, the background ionized fraction in- 
creases extremely rapidly, with most of the ionizations occur- 
ring at z < 4. On the one hand, this means that most bubbles 
form not long before reionization is complete and have rela- 
tively little time to recombine. On the other hand, the fossils 
that do form around the rare luminous quasars at z > 5 have 
an extremely high contrast with the background universe and 
have ample time to recombine nearly fully. 

7.2. A Non-zero Ionizing Background 

We have so far assumed that fossil helium bubbles evolve 
in isolation, without being exposed to any Hell-ionizing ra- 
diation. In contrast to hydrogen reionization, that is not a bad 
assumption - at least for a time - because quasars are so rare, 
and galaxies (probably) do not contribute significantly to he- 
lium reionization. Thus, once a bubble appears, it will remain 
free of any ionizing sources until another quasar forms in the 
same region. The key question is how long this empty phase 
will last. 

To address this question, we must first determine what 
sphere of influence a new quasar will have in the partially 
ionized medium. L^, quasars (with lifetimes of 10^ yr) can 
ionize all the helium inside regions within ^ 15 Mpc of 
their host galaxy. Once regions are mostly ionized, so that 
battling recombinations pr esents the main challeng e (rather 
than ionizing new material). iFurlanetto & OhI (l2007bl) showed 
using the MHROO model that their light is attenuated only 
over ^ 35 Mpc distances. Thus the relevant question is the 
time lag before a quasar appears in this maximum volume 
Knax ~ 30 Mpc. To estimate this, we write the quasar num- 
ber density as nqso = 10^^n_6 Mpc~^ and the lifetime as 
iqso = 10''i7 yrs. For the purposes of a simple estimate, we 
then assume that the quasar population over the time lag At 
can be divided into N discrete generations. The number of 
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Fig. 10. — Ionization liistories inside fossils and in the background uni- 
verse, as in Fig.|2]with Zi = 10. The solid and dashed curves are identical 
to the {zi) = 0.2 and 0.4 curves in Fig.|2]with C = Cmhr- The dotted 
curves also take x^{zi) = 0.2 and 0.4 with C = Cmhr. but they assume 
that photoheating raises the Jeans mass by a factor of ten within the fossils. 



generations required to have a probability of one-half to find 
a second source in the region is then ~ (SnqgoVmax)^^, or 



H{z)At - 0.047' 



tT_ 

n_6 



1 + z 



3/2 



(20) 



According to the iHopkins et alj (l2007h model, the number 
density of quasars with L > 10^^'^ Lq (roughly L^) is n_6 ~ 
(2.2, 0.73, 0.16) at z = (3, 4, 5), implying H{z)At - 
(0.02, 0.06, 0.4:)tj at these redshifts. Thus once reionization 
gets underway, the time lag between successive generations 
of ionizing sources is rather small - only if a region is ionized 
early in reionization will it have sufficient time to recombine 
significantly. For example, a z = 4 (5) bubble will reach 
^Hoiii,™ ~ 0.6 (0.4) before it is reionized by another quasar 
We therefore expect roughly one-half of all photons produced 
at z > 4 to be "wasted" through fossils. But by z = 3.5, the 
fossil lifetimes will be so short that few photons are wasted. 

8. DISCUSSION 
8.1. Photoionization Feedback 

As noted above, it is possible that photoheating increases 
the Jeans mass in the fossil bu bbles (IReesI [19861 lEfstathioul 
ll992tlThoul & Weinberdl 19961) . and with it the effective min- 
imum galaxy mass mmin over that in the universe at large. 
The efficiency of this feedback may be reduced at redshifts 
much larger than z > 6 (Dijkst ra et al.. 2004.) ; nevertheless, it 
can decrease the emissivity in fossils. It is therefore natural to 
ask whether such a reduction may allow fossils to recombine 
more rapidly than we found above. 

Figure [TO] shows a set of ionization histories that address 
this issue. The solid curves are identical to the (z^) — 0.2 
and 0.4 curves in Figure [T] In the dotted curves, we assume 
that photoheating inside the fossil bubble has increased the 
minimum mass to mmin = 1 0m4 (|Rees 1986; Efstathio^ 
119921; iThoul & Weinb erg 1996). Although photoheating de- 
creases the effective collapsed fraction within the bubble by 
nearly an order of magnitude, the effect on xi is relatively 



modest, and the region still remains highly ionized. In part, 
this is because these massive galaxies evolve more rapidly, so 
the emissivity (proportional to d/cou/dz) does not fall by as 
large a factor as J^oW itself. Note also that we still compare 
to a background universe with rrimin — ttla', this is why the 
fossil can have a smaller ionized fraction than the universe at 
large. Of course, in reality photoheating will also slow down 
the background evolution once xf ^, 0.5, an d the two will 
map more smoothly o nto each other (iHaiman & Holder 2003!; 
iFurlanetto & Loeb 2005). But this simple treatment suffices 
to show that photoheating feedback will not allow fossils to 
recombine significantly more than in our fiducial cases; in- 
stead, they are roughly equivalent to a normal history with 
one-half the emissivity. 

However, because massive galaxies are much rarer, inho- 
mogeneity will be easier to re-establish. In particular, Xi 
can be comparable to or smaller than the mean galaxy spac- 
ing if photoheating feedback suppresses galaxy formation in- 
side of fossil HII regions. If we only include galaxies with 
Tvir > 2 X 10^ K (i.e., maximal photoheating feedback), 
dg^i = 4.4, 7.3, 13.7, and 110 Mpc at z 6, 8, 10, and 15. 
Comparing to Figure [TOl the mean free path inside these bub- 
bles can fall below dgai in moderately ionized volumes, even 
at relatively low redshifts. However, the effect of photoheat- 
ing feedback is probably smaller during or just after a region 
is reionized, because objects already in the midst of collapse 
are not very susceptible to it dPiikstra et al.l l2004). Thus it is 
unlikely that suppression at this level can be achieved in the 
relatively short time interval between z = 10 and the end of 
reionization. We would nevertheless generically expect more 
inhomogeneity inside the bubbles if photoheating feedback 
occurs. 



8.2. Detecting Fossil Bubbles 

We have now established that most fossil bubbles formed in 
the midst of hydrogen reionization will remain highly ionized 
throughout the entire process. We will now consider strate- 
gies to find them. To do so, we must rely on two features of 
fossils: their large size, and possibly a different topology of 
ionized gas. We have already seen that the latter effect is not 
so important; most bubbles will remain highly ionized, and 
those that do begin to recombine significantly will resemble 
the background universe but with larger bubbles and less con- 
trast - which will probably only be distinguished with 21 cm 
survey instruments well beyond those currently being built, 
such as the Square Kilometer Array. Instead, it is their large 
sizes that are often taken to be the distinguishing characteris- 
tic of fossils. 

The first potential problem is distinguishing quasar-driven 
fossils from large ionized regions generated by "normal" 
star format ion during reionization. As illustrated in Fig- 
ure |7] ^urlanettoetal] ( I2OO4.) established that HII regions in- 
evit ably reach extremely la rge sizes during reionization (see 
also'Furlan etto et ani2006a ^. During most of the process, the 
characteristic size is relatively small in comparison to the ion- 
ized regions generated by luminous quasars (which are ~ 20- 
40 Mpc across, assuming typical quasar parameters). But this 
is only the characteristic size - in reality, the bubble distri- 
bution has a tail of rare, large objects. Because the quasars 
themselves are also extremely rare, it is not obvious which 
type of rare large bubble will outnumber the other 

Figure [TT] illustrates the problem quantitatively. The three 
thick curves give estimates for the number density of large 
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Fig. 1 1 . — Comparison of the number densities of massive lialos (assumed 
to have hosted luminous quasars) and large bubbles from "normal" reioniza- 
tion (dominated by stellar sources). The thick soHd, long-dashed, and short- 
dashed curves show the abundance of regions with R > 40, 30, and 20 Mpc, 
respectively, from stellar reionization, based on the model of Furlanetto et al. 
(2004). These are evaluated at z = 8 but only depend on redshift slightly. 
The horizontal dotted and solid lines show the total abundance of halos with 
m > 10^'^ M0 and 10^^ Mq, respectively, at the specified redshifts, repre- 
senting possible quasar hosts. 



ionized bubbles from normal star formation.'^ We show 
R = 20, 30, and 40 Mpc, from top to bottom. W e calculate 
these with the analytic model of Furlanetto et alj (.2004): nu- 
merical simulations show that this model underestimates the 
characteristic bubble size early in reionization (by a factor of 
a few), but is quit e accurate in the later stages most applicable 
to this argument (Mesinger & Furlanetto 2007). 

The thin horizontal lines show the number density of halos 
above mass thresholds 10^^ Mq and 10^'^ M fD (solid and dot- 
ted lines, respectively) calculated from the ISheth & TormenI 
mass function at z = 7, 8, and 10 (from top to 
bottom within each set). If we assume that these massive 
halos all host supermassive black holes, then the number 
density of halos above that threshold will equal the maxi- 
mum number density of fossils as well (ignoring recombi- 
nations). The luminous z = 6 quasars are thought to lie 
in halos of m ^ 10^^- 10^'^ Mp., from both number den- 
sity arguments dFan et al.i 2000; Haiman & Loeb 2001) and 
spectral signatures in th e optic al (Barkana & Loeb 20 j|) and 
molecular line spectra (I Walter et al.. ,2004: .Narayanan et al.l 
l2008h . This mass range is also consistent w ith convert- 
ing the empirically-estimated black hole masses (Vestergaard| 
|2004) to a halo mass using the black hole mass-halo mass 
and black hole mas s-velocity dispersion relations calibrated 
at lower redshifts (iMagorrian et al.lll998HFerrarese & MerrittI 
l2000HGebhardt et al.ll2000h ! 

Comparing the bubble and halo number densities, it is clear 
that fossils will not necessarily be distinguishable from more 
normal (galaxy-driven) ionized zones past about the mid- 
way point of reionization. Evidently, the rare tail of ionized 

In detail, we compute Xi(> R)/ (4ttR'^ /3), where Xi{> R) is the 
fraction of space filled by bubbles larger than R. This is more robust to un- 
certainties in the maximum bubble size set by rec ombinations, which bre aks 
extremely large Hll regions into smaller sections IFurlanetto & Ohl2005ft . 



bubbles may provide a comparable number of objects to the 
quasars themselves. That said, the uncertainties in both the 
tail of the bubble distribution (which depends on the overlap 
of smaller regions) and the conversion from quasar luminosity 
to halo mass (and hence number density) are highly uncertain, 
and it is possible that fossils will, in fact, stand out relatively 
well. There are also two morphological effects that could 
be helpful (although difficult to interpret robustly). First, the 
large halos hosting quasars are also surrounded by substantial 
clusters of sources, so they already sit in large ionized bub- 
bles from w hich the quasar c an build a somewhat larger ion- 
ized region (iLidz et al.ll2007h . Second, quasars may present 
distinctive ionization patterns (perhaps two collinear cones, if 
the emission is beamed, or a spherical region if it is isotropic). 

Furthermore, high-redshift galaxy surveys over the same 
fields may allow us to break this degeneracy. 21 cm 
and galaxy surveys should have a strong detectable anti - 
correlation dFurlanetto & Lidzl 120071: IWvithe & Lo"ebll2006h . 
Moreover, on the largest scales where bubbles are directly de- 
tectable, we could simply compare the observed comoving 
emissivity within the bubble to its size.'* This relation will 
have some natural scatter due to the past star formation his- 
tory within the bubble, variations in gas clumping, and ob- 
servational error However, since quasars constitute a strong 
impulsive spike in the comoving emissivity, fossil quasar bub- 
bles will induce a decorrelation between the observed comov- 
ing emissivity and bubble size. In particular, quasar fossils 
will lie far off the tail of the emissivity distribution for a given 
bubble size (with the deviation being proportional to the frac- 
tion of ionized volume contributed by the quasar). Such ef- 
fects could potentially allow us to distinguish quasar-blown 
bubbles (the majority of which will be fossil regions), permit- 
ting a lower bound on the contribution of quasars to reioniza- 
tion. 

Another potential problem in identifying these bubbles is 
that the underlying ionization (and density) fields are fluctu- 
ating on large scales from the cumulative effects of smaller 
regions. We need to be able to distinguish the fossil region's 
signal from this background "confusion noise". The best way 
to discover fossil bubbles will probably be through the 21 cm 
line dFurlanetto et al. 2006b), for which the brightness is pro- 
portional to the total density of neutral gas. A fully ionized 
fossil would have a contrast ^ 20a;Hi mK with the mean back- 
ground. To estimate the conf usion limit in such an observa- 
tion, we use the simulations of Mesinger & Furlanettol (l2007h 
to measure the r.m.s. fluctuation amplitude at wavenumbers 
k ^ (27r/30) Mpc^^. This varies from ~ 0.1-1 mK over 
the range xhi ~ 0.8-0.1. Thus we would expect fossils to 
stand out very well - many standard deviations from the mean 
- except near the end of reionization, when they will begin 
to blend into the background fluctuations as well. The main 
obstacle toward direct detection is therefore telescope noise, 
though sufficiently large bubbles may still be detectable, even 
by the first generation of instrumen ts, at least at low redshifts 
(see §8.6 of lFurlanetto et al.ll2006bl) . 

Although the uncertainties in these estimates are fairly large 
(especially in the abundance of massive quasars), these fac- 
tors do suggest that it may be significantly more difficult to 
cleanly identify fossil bubbles than often assumed, at least 

Note that the bubble size is proportional to the integrated star formation 
history oc /coii , rather than the instantaneous comoving emissivity oc /coii ■ 
However, the two are strongly correlated in the sense that the rarest high den- 
sity regions with the highest /coii will also have the highest rate of collapse. 
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late in reionization: they may not be much more abundant 
than large bubbles from star formation, they may not stand out 
strongly from the background field, and their ionization pat- 
tern (either remaining highly ionized or with discrete ionized 
regions separated by partially ionized walls) does not qual- 
itatively differ from the background. We suspect that some 
modeling will be required to distinguish the two cases. 

8.3. Fossil and Active Bubbles 

As we have seen above, although the fossils can be kept 
highly ionized, they are not fully ionized, and will thus ap- 
pear different from active quasar bubbles (for example, in 2 1 
cm imaging, they will have a somewhat reduced contrast). 
Since active quasars, with the luminosities required to pro- 
duce large (> 10 Mpc) fossil bubbles, should be relatively 
easily detectable even at z > 6, active bubbles can probably 
be distinguished routinely from fossil bubbles, simply by the 
presence or absence of a luminous ionizing source. 

We now turn to the question of how many fossil bubbles 
we might expect relative to those surrounding active quasars. 
Fossils forming relatively close to the end of reionization (at 
z < 10 for reionization at z 7) do not significantly re- 
combine so would appear as mostly empty fossils long after 
their source quasar fades. This will boost the number of large 
(mostly) ionized regions well above that expected from the 
quasar luminosity function alone, an d can potentially allow 
one to constrain th e quasar lifetime jWvithe & Loebll2004bl : 
IWvithe et al.ll20Q5h . 

However, because fossils persist for so long, the quasar 
abundance itself probably evolves significantly over the re- 
combination time. As an example, let us suppose that the 
number density of bright quasars follows nqso oc (1 + z)^^^ , 
with /3 > 0, and that fossils remain ionized so long as they 
form at z < z,„ax- In that case, assuming that the quasar life- 
time iqso is much smaller than the elapsed time from Zmax 
to the observed time Zobs, we expect 

'^QSo(^obs) 



^fossil 
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{13 + 3/2)i?(zob«)tQso ■ 
The number density of bright quasars has /3 ^ 4 over the 
range z ^ 3-6 (Fanetal. 2001, 2004). If this continues to 
higher redshifts, it will provide a factor ^ 5 suppression from 
the naive value. Unless it can be measured independently, 
this evolutionary factor will pose a significant (factor of sev- 
eral) uncertainty in the measurement of iqso from the fos- 
sil to active ratio. One way to break the degeneracy is by 
measuring the luminosity function over the entire time inter- 
val - but doing so would require a search for active quasars 
that extends to much higher redshifts than the fossil search 
itself. Another potential method would be a series of mea- 
surements of the number density of bubbles over a redshift 
interval comparable to the recombination time of a bubble. 
This could probe the redshift evolution of the quasar lumi- 
nosity function, particularly if sensible priors for the quasar 
lifetime are adopted. Finally, we note that the shapes of ac- 
tive quasar bubbles could, in principle, contain information on 
the hos t quasars and the IGM, through finite light-travel time 
effects (IWyithe & Loeb 2004b; Yu 2005). The results we ob- 
tained here, i.e. that most of the fossil gas takes significantly 
longer to recombine than one would estimate using a fidu- 
cial IGM clumping f actor, will ham per statistical versions of 
this measurement ( Sethi & Haiman 2008), since fossils will 
outnumber active quasar bubbles and not show the apparent 
distortion due to finite light travel time. 



8.4. Double Reionization 

Up to this point, we have studied recombinations in the con- 
text of fossil bubbles, but much of our discussion also applies 
to so-called "double r eionization" scenarios CCen 2003bja|; 
IWvithe & Loebll2b03allbh in which a first generation of mas- 
sive, metal free stars reionizes most of the universe but also 
enriches it - inducing a transformation to much less efficient, 
but more normal, stars. In the scenarios originally proposed, 
photons from the second generation did not suffice to bat- 
tle recombinations and so the universe became substantially 
(up to 80%) neutral, until structure formation had progressed 
enough for the more normal stars to complete reionization. 

However, these models rely on relatively large clumping 
factors to speed up recombinations so that a recombination 
era can occur during the relatively short interval between the 
two ionizing source generations. Reionization requires that 
(at least) one ionizing photon per hydrogen atom be produced; 
this is relatively difficult at high redshifts. But maintaining a 
high ionization fraction is less stringent: it requires only that 
one ionizing photon be produced per recombination. For gas 
at the mean density, where the recombination time is roughly 
the Hubble time, this is relatively easy. We have seen that the 
effective clumping factor declines as more gas recombines, 
suggesting that the "recombination era" will probably be con- 
fined to dense pockets of gas. 

Our models strengthen this conclusion. Figure |7] can be 
used to analyze double reionization scenarios in addition to 
fossil bubbles: we can simply think of the if (z) curves as 
describing the second generation, and we suppose that the 
fossil's initialization redshift describes the moment at which 
the first generation shuts off. The left and right panels then 
mimic moderate and more extreme double reionization sce- 
narios. Our most important conclusion is that, even in the 
most extreme case - when the second generation is able to 
achieve only Xi ^ 0.005 at the initial point - less than half 
of the gas is actually able to recombine. In more moderate 
scenarios, only ~ 10% can recombine, and that is generally 
dense gas surrounding galaxies and halos. 

Note finally that these scenarios assume a sudden tran- 
sition between the two generations - in reality, that is 
much too simplistic, because all of the feedback mech- 
anisms that can mediate such a transition (metal enrich- 
ment, photoheating, etc.) act locally and on relatively long 
time scales (Haiman & Holder 2003; Diikstra et al. 20()4 
Furla netto & Loebl 120051; llliev et al.1 |2007; Mc Ouinn et al] 
i2007i) . Both of these effects (slow recombinations and slow 
transitions) will combine to make double reionization excep- 
tionally unlikely. 

9. SUMMARY 

We have examined recombinations in inhomogeneous fos- 
sil ionized regions, adapting the simple MHROO model of the 
density distribution of the IGM, and tracking the ionization 
histories of individual gas elements. Without a residual ion- 
izing background, the process occurs quickly at first, because 
the dense IGM gas can recombine rapidly (during both hy- 
drogen and helium reionization). But once these relatively 
rare pockets become neutral, recombinations slow dramati- 
cally, and eventually Ccfr < 1. Thus it is quite difficult for a 
region to recombine entirely, even in the absence of any ion- 
izing flux. 

The ionizing background from galaxies inside the fossil, to- 
gether with any residual low-level emission from the black 
hole past its bright quasar phase, efficiently further suppresses 
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recombinations during hydrogen reionization. So long as 
galaxies are able to ionize > 20% of the universe on average, 
they will be able to halt recombinations in all but the dens- 
est gas inside the fossil regions. In most cases, the mean free 
path of ionizing photons remains large enough to wash out any 
fluctuations in the ionizing source population, and the fossil 
remains uniformly ionized (outside of self-shielded regions). 
In more extreme scenarios, when the fossil is created long be- 
fore reionization, the gas can recombine significantly (down 
to roughly the mean density), which wiU probably allow large 
regions far from galaxies to shield themselves and begin re- 
combining faster. In such a case, the fossil IGM will look 
quahtatively similar to the swiss-cheese ionization topology 
of the IGM elsewhere, but with larger ionized regions and a 
smaller contrast between the ionized and "gray" partially neu- 
tral regions. 

These fossil regions may therefore be somewhat harder to 
distinguish from the normal IGM than often assumed - more 
so because "normal" ionized regions can reach comparable 
sizes to those of quasar zones near the middle of reionization, 
and may be as or even more abundant (albeit both of these 
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estimates rely on the rare tails of theoretical distributions, in 
galaxies and ionized bubbles, and have very large uncertain- 
ties). 

During helium reionization, fossil bubbles evolve some- 
what differently. Because the ionizing sources are rare, once 
a region's source quasar shuts down, it is unlikely to be illu- 
minated by any other source for a time. This will allow the 
gas to recombine without interference; because the universe 
is so clumpy at z ~ 4, a large fraction of the gas can recom- 
bine over that interval. However, the typical time lag before 
the next quasar appears is only a fraction of the Hubble time, 
so in practice less than half of the gas will be able to recom- 
bine, even for quasar bubbles formed early in the reionization 
process. 
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